Spin Hall Effect 



John Schliemann 
Institute for Theoretical Physics, University of Regensburg, 
D-93040 Regensburg, Germany 

February 4, 2008 



Abstract 

The intrinsic spin Hall effect in semiconductors has developed to a remarkably lively and 
rapidly growing branch of research in the field of semiconductor spintronics. In this article 
we give a pedagogical overview on both theoretical and experimental accomplishments and 
challenges. Emphasis is put on the the description of the intrinsic mechanisms of spin Hall 
transport in III-V zinc-blende semiconductors, and on the effects of dissipation. 

1 Introduction 

Starting from about the late 1990s, a pronounced, and partially also renewed, interest in effects 
of spin-orbit coupling in semiconductors has been emerging in the solid-state community. This 
development is mainly fueled by the field of spintronics. The latter keyword summarizes an 
entire plethora of theoretical and experimental efforts towards using the spin degree of freedom 
of electrons, instead, or in combination with, their charge for information processing, or, even 
more ambitious, for quantum information processes. Thus, controlling the electron spin in 
semiconductor structures is a key challenge, and the relativistic effect of spin-orbit coupling 
is an important, if not indispensable, ingredient for reaching this goal. Among all the rapidly 
progressing activities in this field, a major and very recent development was the theoretical 
prediction and subsequent experimental investigation of the spin Hall effect in semiconductor 
structures. This effect amounts in a spin current, as opposed to a charge current, driven by a 
perpendicular electric field. In this article we review recent studies on spin Hall transport in 
semiconductors induced by intrinsic mechanisms. 

A brief overview on important aspects and perspectives of the field of spintronics is given 
in the article by Wolf et al. [1] . Selected topics are reviewed in more detail in a volume edited 
by Awschalom, Loss, and Samarth [2]. A comprehensive review of many parts of spintronics 
was given by Zutic, Fabian, and Das Sarma [3]. Research directions not covered here include 
the field of ferromagnetic semiconductors; for a review on this we refer to Refs. [4, 5, 6, 7, 8, 9]. 

From a historical perspective, the notion of the spin Hall effect in systems of itinerant 
spinful charge carriers was considered first by Dyakonov and Perel [10] in the early seventies, 
and in a more recent paper by Hirsch [11]. In these studies the predicted spin Hall effect is 
due to spin-orbit effects influencing scattering processes upon static impurities. Following the 
usual terminology of semiconductor physics, this effect is referred to as the extrinsic spin Hall 
effect since it necessarily requires spin-dependent impurity scattering. This is in contrast to 
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the intrinsic spin Hall effect which is entirely due to spin-orbit coupling terms in the single- 
particle carrier Hamiltonian and occurs even in the absence of any scattering process. The 
story of the intrinsic spin Hall effect starts in summer of 2003, when Murakami, Nagaosa, 
and Zhang [12], and, almost simultaneously, J. Sinova et al. in a collaboration based in 
Austin (Texas) [13] predicted this phenomenon. The work by Sinova et al. considers a two- 
dimensional electron gas being subject to spin-orbit coupling of the so-called Rashba type, 
whereas the former paper investigates valence-band holes in three-dimensional bulk systems. 
Shortly later on, Schliemann and Loss added a study on spin Hall transport of heavy holes 
confined to a quantum well[14]. Only a few months later, J. Wunderlich et al. reported on 
an experimental study of the latter type of system, confirming the predicted spin Hall effect 
via optical techniques [15]. The paper by Wunderlich et al. was only shortly preceded by 
another experimental report by Kato et al. who detected spin Hall transport in n-doped bulk 
systems, again using an optical method [16]. In both experiments, the existence of spin Hall 
transport is signaled by spin accumulation at the boundaries of the sample, which seems to 
be the easiest method so far to detect this effect. 

The theoretical and experimental developments sketched above have generated a still 
rapidly growing amount of preprints and (subsequent) journal publications, so far mostly 
theoretical. In this article we give a pedagogical overview on important aspects of intrinsic 
spin-Hall transport in III-V zinc-blende semiconductors. Emphasis is put on the the descrip- 
tion of the intrinsic mechanisms of spin Hall transport induced by spin-orbit coupling,and 
on the effects of dissipation. This article also addresses researchers who are not particularly 
specialized in spin phenomena in semiconductors but are interested in this rapidly develop- 
ing field. A brief review on spin-Hall transport was already given in a conference paper by 
Murakami [17], and Sinova et al. recently provided a short summary of important issues [18]. 

This article is organized as follows. In section 2 we make a few general comments on spin- 
orbit coupling in semiconductors and describe its effective contributions to the band structure 
of both electron and hole doped systems. In section 3 we first make some general remarks 
on the notion of spin currents before reviewing the particularly rich body of recent theoret- 
ical work on spin Hall transport. Our analysis here includes the two-dimensional electron 
gas as well as p-dopcd bulk systems and quantum wells. Experimental work and proposed 
experiments are discussed in section 4. We close with conclusions and an outlook in section 5. 

2 Spin-orbit coupling in III-V semiconductors 

The coupling between the orbital and the spin degree of freedom of electrons is a relativistic 
effect described by the Dirac equation and its nonrelativistic expansion in powers of the inverse 
speed of light c. In second order one obtains, apart from two spin-independent contributions, 
the following well-known spin-orbit coupling term. 



where mo is the bare mass of the electron, s, p its spin and momentum, respectively, and V is 
some applied external potential. On the other hand, the free Dirac equation, V = 0, has two 
dispersion branches with positive and negative energy. 
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which are separated by an energy gap of 2mQC^ ~ IMeV. In particular, the nonrclativistic 
expansion of the Dirac equation quoted above can be seen as a method of systematically in- 
cluding the effects of the negative-energy solutions on the states of positive energy starting 
from their nonrelativistic hmit. Moreover, the large energy gap 2moc^ appears in the denomi- 
nator of the right hand side of Eq. (1), suppressing the effects of spin-orbit coupling for weakly 
bound electrons. 

On the other hand, the band structure of zinc-blende III-V semiconductors shows many 
formal similarities to the situation of free relativistic electrons, while the relevant energy scales 
are grossly different [19, 20, 21]. For not too large doping of such semiconductors, one can 
concentrate on the band structure around the T point. Here one has a parabolic s-type 
conduction band and a p-type valence band consisting of the well-known dispersion branches 
for heavy and light holes, and the split-off band. However, the gap between conduction and 
valence band is of order leV or smaller. This heuristic argument makes plausible that spin- 
orbit coupling is an important effect in III-V semiconductors which actually lies at the very 
heart of the field of semiconductor spintronics. 

In the following we give an overview on effective model Hamiltonians for conduction- 
band electrons and valence-band holes in III-V zinc-blende semiconductors in several spatial 
dimensions. These effective expressions can be obtained via the so-called k ■ p-theory and 
related methods, as general references we refer to Refs. [22, 23, 24]. Here we shall just state 
the results and discuss their main physical implications. 

2.1 Conduction-band electrons 

Let us first consider three-dimensional bulk systems. For electrons in the s-type conduction 
band, the contribution to spin-orbit coupling being of lowest order in the electron momentum 
p has been derived by Dresselhaus [25] and reads 

where a is the vector of Pauli matrices describing the electron spin, and 7 is an effective 
coupling parameter. This Hamiltonian is trilinear in the momentum p and invariant under 
all symmetry operations of the tetrahedral group T^, the point symmetry group of the zinc- 
blende lattice. As a result, the parameter 7 is different from zero because the zinc-blende 
lattice does not possess an inversion center. Therefore, the Dresselhaus spin-orbit coupling is 
due to bulk-inversion asymmetry. 

In a sufficiently narrow quantum well grown along the [001] direction, and at sufficiently 
low temperatures, one can approximate the operators Pz and p^ by their expectation values 
[pz) « 0, (pI) = Ti^{k^). Then neglecting terms of order p^, Py leads to a spin-orbit coupling 
term linear in the momentum [26, 27], 

nD = ^{Pyay-p,a-) (4) 

with P = ^{k1). Here kz is the wave number in the lowest subband of the well. Another 
important contribution to spin-orbit coupling occurs in quantum wells whose confining poten- 
tial is lacking inversion symmetry. This contribution due to structure-inversion asymmetry is 
known as the Rashba term [28, 29], 

nR = ^{p.ay-pyan, (5) 
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where the coupUng parameter a is essentially proportional to the potential gradient across 
the quantum well. Hence, a is in particular tunable by an electric gate and can therefore be 
varied experimentally. 

A both theoretically [30] and experimentally [31, 32] well established value for the Dres- 
selhaus parameter in GaAs is 7 = 25eVA^. Depending on the width of the quantum well, this 
leads to values for [5 being of up to 10~^^eVm. Regarding the Rashba coefficient a, values 
of a few 10~^^cVm can be reached in InAs [33, 34, 35, 36, 37, 38, 39], whereas in GaAs this 
quantity is typically an order of magnitude smaller [40]. Thus, the characteristic energy scales 

£R = -J^ > (6) 



mfP 
1? 

for Rashba and Dresselhaus coupling, respectively, can be of order 0.1 .. . l.OmeV, depending 
on the effective band mass m. 

Let us finally briefly discuss the spectrum and eigenstates generated by the above spin-orbit 
coupling terms. We consider the single-particle Hamiltonian for a two-dimensional electron 
system 

n = ^ + nR + nD. (8) 

The eigenenergies are given by 

e± {k) = '^± ^Jiaky + f3k^f + (afc^ + f5kyf (9) 
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with eigenstates 

where A is the area of the system and 

X{k) = arg{-aky - f3kx + i{akx + (3ky)) . (11) 

The above spin-orbit coupling terms can be viewed as a momentum-dependent Zeeman field 
acting on the electron spin. Consequently, the spin state of the electron depends on its 
momentum, as seen in Eq. (10). Note that for pure Rashba or Dresselhaus coupling, the 
dispersions form two parabolas being shifted horizontally. This is different form a normal 
Zeeman field which shifts the dispersion parabolas vertically, i.e. along the energy axis. The 
case a = ^(3 is particular [41, 42]. Here a new conserved quantity given byS := ((T^TCT^')/^ 
arises, and the spin state of the electrons becomes independent of the wave vector. This result 
for a = ib/3 is a very general one, it also holds in the presence of any arbitrary scalar potential, 
or if interactions between electrons are included. 

2.2 Valence-band holes 

The valence band of HI-V zinc-blende semiconductors is of p-type, i.e. it is predominantly 
composed out of atomic wave functions with angular momentum 1 = 1. Adding this angular 
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momentum to the electron spin s = 1/2, we find to multiplets with total angular momentum 
j = 3/2 and J = 1/2. The dublctt j = 1/2 forms essentially an energetically separated so-called 
split-off band and will not be considered any further. The multiplet j = 3/2 consists essentially 
of the so-called heavy and light hole states which are, to a good degree of approximation, 
described by Luttinger's Hamiltonian [43], 

Here iuq is again the bare electron mass, and S are spin-3/2-operators. The dimensionless 
Luttinger parameter 71 and 72 describe the valence band of the specific material with effects 
of spin-orbit coupling being included in 72. The eigenstates of the above Hamiltonian can be 
chosen to be eigenstates of the helicity operator X = (k ■ S)/k, where k = p/h is the hole wave 
vector. The heavy holes correspond to A = ±3/2, while the light holes have A = ±1/2. Prom 
the Hamiltonian (12) one finds the effective band mass for the heavy holes as 



and for the light holes 



ruhh = 7^ (13) 

71 - 272 

mh = — —TT- ■ (i^^) 

71 + 272 

Well established values for the Luttinger parameters, among other band structure parameters, 
can be found in the literature [44]. For example, for GaAs one has 71 ~ 7.0 and 72 ~ 2.5 
giving nihh ~ 0.5mo and mih ~ O.OSmo- 

In a bulk system, the heavy and light hole states are degenerate at the F-point k = 0. This 
degeneracy is lifted in a quantum well due to size quantization, and for sufficiently narrow wells 
and low enough temperatures one can concentrate on the lower-lying heavy holes. Moreover, 
if the quantized wave vector in the growth direction is large enough, i.e. if the well is not 
too wide, the spin of these heavy holes points predominantly along the growth direction 
with a projection of ±3/2. For asymmetric wells these hole states are subject to a spin-orbit 
contribution due to structure-inversion asymmetry analogous to the Rashba term for electrons 
in the conduction band. Choosing the growth direction to point along the 2;-axis, the resulting 
effective Hamiltonian has the form [45, 46] 

using the notations p± = Px ^ iPy, cr± = i: ia^. Here the Pauli matrices operate on the 
total angular momentum states with spin projection ±3/2 along the growth direction; in this 
sense they represent a pseudospin degree of freedom rather than a genuine spin 1/2. In the 
above equation, m is the heavy-hole mass, and a is Rashba spin-orbit coupling coefficient due 
to structure inversion asymmetry. This Hamiltonian has two dispersion branches given by 

e±{k) = -^±ak^ (16) 



with eigenfunctions 
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We note that the Rashba parameter a entering the Hamiltonian (5) for electrons in a quantum 
well has a different dimension than the parameter a for holes in Eq. (15). In the latter case, 
the quantity majf? has dimension of length and can reach a magnitude of several nanometers 
in GaAs samples [47]. 



3 Spin Hall transport: Theory 

In this section we summarize important theoretical results on intrinsic spin Hall effect. We 
start with some general considerations on spin currents. 



3.1 Spin currents: General remarks 

A type of current most familiar to physicists is certainly the usual charge or particle current. 
The particle density in a many-body system is described by the operator 

p{?)=Y.^{f-?n) (18) 

n 

where the index n labels the particles. This density operator is a function of time via the 
time-dependence of the positions r„ entering the argument of the delta-functions. Let this 
time evolution be generated by an Hamiltonian of the form 

'H = ^h{fn,Vn-,Sn)^yint- (19) 
n 

Here the term Vint describes interaction among the particles and depends only on their spatial 
coordinates, and the single-particle Hamiltonian h reads 

^P) = ^ + lijVi^' + V^(r) . (20) 

Summation over repeated cartesian indices is understood, and the matrix 7 parameterizes spin- 
orbit coupling of the Rashba and Dresselhaus type for electrons in a quantum well. Finally, 
the static potential Vif) describes, e.g., static impurities. Now, starting from the Heisenberg 
equation of motion, 

and performing some elementary algebraic manipulations, one derives the well-known conti- 
nuity equation for the particle current, 

|p + V.J=0 (22) 

with the current density operator 

J(0 = \ '^{^{Vn,dn),^{f - rn)} ■ (23) 
^ n 

Here {^, B) = AB + BA denotes the anticommutator of two operators, and the velocity 
operator v is given by 

v{Pn,^n) = ^[H,fn] (24) 
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for each particle n. Note that this operator is in general spin-dependent if spin-orbit coupling 
is present. 

Let us now consider a general observable described by a hermitian single-particle operator 
A which can be a function of position, momentum, and spin: A = A{f,p,a). The density 
operator corresponding to this physical quantity A is naturally defined as 

PA{r) = ^'^{A{'rn,Pn,^n),^{'r-fn)} , (25) 

n 

where the symmetrization ensures hermiticity. Now proceeding as above one finds 

^PA + V • = , (26) 
where the current density operator Ja is given by 

jA{r) = ^^{A{rn,Pn,ffn),{v{Pn,^n),Sir-fn)}} , (27) 
n 

and the additional source term on the right-hand-side reads 

SAifO = lY.{^i^^'^(''n,Pn,ffn)],Sif-fn)^ . (28) 

Thus, wc only obtain the usual from of the continuity equation if the observable A commutes 
with the Hamiltonian Ti., which is of course just a restatement of Noether's theorem. For 
the case electron spin components as observables, the corresponding spin-current densities are 
given by 

U^ = ]E{^<dviPn,an),S{f-fn)}\ , (29) 
n ^ ' 

and the source terms are due to spin-orbit coupling being present in the single-particle Hamil- 
tonian. These source terms reflect the fact that magnetization, i.e. the density of magnetic 
moments, can be altered by two ways: by spatially moving spinful particles, or by manipulat- 
ing their spin state. The latter process is described by the source terms. For instance, in the 
case of noninteracting electrons in a quantum well with spin-orbit coupling of the Rashba and 
Dresselhaus type, the source terms can be expressed via the components of the spin-current 
densities itself [48, 49], 



d „ 2ma „ 2mf3 ., 



fc+V-j, = — jf + -fj», (30) 



^P. + V.,, = -^,t + -^jt, (31) 

|p.+V.i = (32) 

In summary, the definition of the spin current density as given in Eq. (29) is the straightforward 
generalization of the usual particle current and widely used in the literature. As seen above, 
this spin current density is, however, not conserved, i.e. it does not fulfill a simple continuity 
equation. This fact might or might not be seen as a shortcoming of the above definition. 
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Another peculiarity of this type of current operator was pointed out by Rashba [50] who 
found that, for the situation of an asymmetric quantum well, the current densities with in- 
plane spin components, jx, jy, can have nonzero expectation values even in the absence of an 
electric field, i.e. in thermal equilibrium. Using periodic boundary conditions and considering 
an infinite disorder-free system of non-interacting electrons at zero temperature and positive 
Fermi energy, the full result for the case of both Rashba and Dresselhaus coupling reads [48] 

Note that these equilibrium spin currents vanish in the case a = ib/3 due to the additional 
conserved spin operator arising at this point [41]. The findings shown in Eqs. (33), (34), 
however, certainly depend on the boundary conditions used and are altered in a more realistic 
description of finite systems [51]. 

In the recent literature, there are several proposals and discussions on alternative forms of 
spin currents which possibly fulfill proper continuity equations [52, 53, 54, 55, 56, 57]. However, 
these issues do not seem to be settled yet. Therefore, in the following we shall concentrate on 
spin current densities as defined in Eq. (29). 

3.2 Conduction-band electrons in two dimensions 

We now discuss spin Hall transport of conduction-band electrons in III-V semiconductor quan- 
tum wells. As a great simplification used in almost the entire theoretical work so far, we will 
consider non-interacting electrons. Thus, the system is described by the single-particle Hamil- 
tonian 

W = ^ + f iP.^' - Py<Tn + f {Py^y - pxa-) , (35) 

and instead of the many-body spin current density operators (29) we can use the single-particle 
operator 

Jz = j{a'v + va') (36) 
m 2 

where we have concentrated on the spin component along the growth direction of the quantum 
well (chosen as z-axis) and used the anticommutativity of Pauli matrices. To account for effects 
of disorder and confining boundaries of the system, appropriate potentials should be added to 
the Hamiltonian (35), as we will discuss in detail below. 

The linear response of this spin current to an electric field applied in the plane of the 
two-dimensional electron gas can be evaluated via the usual Kubo formula [58]. For the 
off-diagonal (or Hall) components of the response tensor one has [13, 59] 



p poo 

('"^ = A{J+ir^) Jo «^^"^^''^*E/(^MW)(^>)"lb-fW,^.(0)]|/c,/x). (38) 
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Here A is the volume of the system, e is the elementary charge, and f{e^{k)) is the Fermi 
distribution function for energy e^{k) at wave vector k in the dispersion branch ^ = it as 
given in Eq. (10). The above quantity describes the linear response in terms of a spin current 
to a perpendicular electric field of frequency 00. In the commutator on the right-hand side the 
time-dependent spin current operator in the Heisenberg picture enters, 

J,(^) = e^^*/'^J,e-^«*/^ (39) 

Moreover, tlic right-hand side of Eq. (38) lias to be understood in the limit of vanishing 
imaginary part 77 > in the frequency argument. This imaginary part in the frequency reflects 
the fact that the external electric field is assumed to be switched on adiabatically starting from 
the infinite past of the system, and it also ensures causality properties of the retarded Green's 
function occurring in Eq. (38). In general, and as we will discuss in more detail below, the 
limiting process 77 — > does not commute with other limits, and, in particular, the dc-limit 
a; — has to be taken with care [58]. In the presence of random impurity scattering, the 
retarded two-body Green's function in Eq. (38) will generically have a frequency argument 
with positive imaginary part [58]. In this case the limit r/ ^ is unproblematic, and the 
imaginary part of the frequency argument is just due to impurity scattering and/or other 
(many-body) effects. Generically, the imaginary part ry > corresponds to a finite carrier 
quasiparticle lifetime. 

Let us now for simplicity consider the case of Rashba spin-orbit coupling only and assume 
the electron density to be large enough such that the Fermi energy is positive (which is usually 
the case in realistic samples) . Neglecting all possible disorder effects and concentrating on the 
case of zero temperature, the Kubo formula (38) can be evaluated straightforwardly, giving a 
spin Hall conductivity at zero frequency of 

^fr(0) = -^^'(0) = ^- (40) 

This result was obtained first by Sinova et al. [13]. It is remarkable in the sense that the value 
of the spin Hall conductivity does not depend on the Rashba parameter a. In particular, 
even in the limit of vanishing spin-orbit coupling, the above result still predicts a finite spin 
Hall current. However, no effects of disorder in the system have been included so far, and 
the infinitesimal parameter rj has been put to zero right away. Let us now take into account 
disorder effects by replacing 77 with a phenomenological relaxation rate 1/r. Here we find [59] 

.f.(0).-.-,0).^-^^...-(4i^(l + 8^)-) . ,41) 

The first term is still the universal expression found in Ref. [13], whereas in the second contri- 
bution three energy scales enter: The Fermi energy ej, the Rashba energy e/j, and the energy 
scale of the scattering by disorder potentials K/t. Clearly, if the latter quantity dominates 
over the Rashba coupling, ?i/T » er, the second term in Eq. (41) cancels the first one, and 
the spin Hall conductivity is indeed suppressed by disorder. Analogous results can be found 
if both the Rashba and the Dresselhaus coupling are included [60]. In this case, the above 
approach also yields nonvanishing longitudinal spin conductivities [60] . 

Thus, we arrive at an apparently physically satisfactory picture. However, it turns out to 
be qualitatively incorrect for the following reason: When replacing the infinitesimal parameter 
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rj in Eq. (38) with the inverse of a phenomenological relaxation time r, one neglects certain 
contributions in a systematic perturbational expansion with respect to the disorder potentials. 
These contributions are known as vertex corrections. As is was shown first by Inoue, Bauer, 
and Molenkamp [61], the full dissipative contribution to the spin-Hall conductivity including 
the vertex corrections exactly cancels the universal value, independently of the strength of 
the disorder potentials and the Rashba coupling. This result was obtained for an infinitely 
large system and in lowest perturbational order with respect to the disorder potentials which 
were modeled by delta-functions. Subsequently, this finding was reproduced an generalized 
by several other authors [62, 63, 64, 65, 66, 67, 68] using different theoretical methods. The 
conclusions from these investigations can be summarized as follows: The spin-Hall conductivity 
for spin polarization along the growth direction in an infinite two-dimensional system with 
spin-orbit coupling of the Rashba and Dresselhaus type strictly vanishes in the presence of any 
spin-independent mechanism with forces, via spin-orbit coupling, the electron spins to relax 
to a constant value. This result is independent of any perturbational expansion with respect 
to disorder terms and holds also both at finite temperature and in the presence of interactions 
among the electrons. However, it does in general not hold in the presence of magnetic fields 
or other spin-dependent contributions in the Hamiltonian. 

The proof of this very general statement was worked out by Chalaev and Loss [65], and 
by Dimitrova [66]; a preliminary version can also be found in Ref. [48]. For the case of 
both Rashba and Dresselhaus spin-orbit coupling, the argument is as follows: To analyze the 
electron spin dynamics in an infinite homogeneous system, it is sufficient to consider just the 
spin operator of a single electron. Prom the Heisenberg equation of motion we find 



d „ 4ma „ AmjS ... , 

-rcr = — H T-fz , (42) 

dt Jz ^3 Jz, \ J 



These relations hold also in the presence of any arbitrary spin-independent potential or in- 
teraction term in the Hamiltonian. The key observation is that the time derivatives of the 
spin components can be expressed as linear combinations of the spin current operators itself. 
This result crucially relies on the fact that the spin-orbit coupling is linear in the electron 
momentum. Moreover, in the presence of spin-orbit coupling disorder effects can generally be 
expected to make the electron spins relax to constant values. Thus, in a stationary state, the 
expectation values of the left-hand sides of Eqs. (42), (43) should vanish. Now it follows im- 
mediately that the expectation values of the spin current components , must also vanish 
provided a / ±/?. In the case a = ±/3, however, spin Hall transport is generally absent due 
to the additional conserved spin quantity which occurs here [60, 69, 41]. 

The central argument of the above proof is closely related to studies by Erlingsson, Schlie- 
mann and Loss [48], and by Dimitrova [66] on the relationship between spin currents and 
magnetic susceptibilities. For further developments in this direction see also Ref. [70]. 

Another interesting observation regarding spin Hall transport in n-doped quantum wells 
was made by Rashba [71] who considered Rashba spin-orbit coupling in the presence of a 
magnetic field coupling to the orbital degrees of freedom only, neglecting the Zeeman coupling 
to the electron spin. In this model, the spin Hall conductivity vanishes in the limit of vanishing 
magnetic field even in the absence of disorder, an effect closely related to the abovementioned 
vertex corrections [71]. Thus, the case of zero magnetic field (coupling to the orbital degrees 
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of freedom only) is different from the limit of vanishing field. This result certainly relies on 
the assumption of an infinite system, since the orbital effects of a magnetic field should be 
small if the typical cyclotron radius is large compared to the system size. 

The case of a magnetic field coupling both to the orbital and spin degree of freedom of 
electrons was investigated by Shen et al [72, 73]. Due to the coupling to the spin this situation 
is not covered by the above general argument. Shen et al find a resonant behavior of the 
spin Hall conductance as a function of the magnetic field when a degeneracy of Landau levels 
occurs at Fermi level. These studies, however, do not include disorder effects so far [72, 73]. 

Coming back the case zero external magnetic field, Adagideli and Bauer chose yet another 
approach to the problem by considering the electron acceleration in the presence of Rashba 
coupling [74], 

^r = ^,.x4--(eE + VF), (44) 

where is the unit vector in the z-direction, E is the in-plane electric field, and V is the 
disorder potential. Since the acceleration should relax to zero in a disordered system, the 
spin Hall current vanishes if the effects of the electric field and the disorder potential cancel 
on average. As the authors show, this is indeed the case in the bulk of the system, but not 
necessarily at its edges. This observation gives rise to the notion of spin Hall edges [74]. 

The above general result on the absence of spin Hall transport in infinite systems with spin- 
orbit coupling being linear in the electron momentum is the outcome of an intense theoretical 
discussion in the last about two years. It was also confirmed by numerical studies carried out 
by Nomura et al. [75]. These authors performed a careful numerical evaluation of the Kubo 
formula for the spin Hall conductivity in the presence of Rashba coupling and delta-function 
type impurity potentials. 

We stress again that the above general conclusion rules out spin Hall transport only in 
the limit of an infinite system. In fact, several numerical investigations on finite systems 
have appeared recently [76, 77, 78, 79, 80, 81, 82, 83, 84]. In these studies, the underlying 
semiconductor structure are described by tight-binding hopping models of finite-size coupled 
to semi-infinite leads. These hopping Hamiltonians also include local disorder potentials and 
discrete versions of the spin-orbit contributions in an n-doped quantum well. Transport quan- 
tities are then evaluated using the well-established Landauer-Biittiker approach combined with 
a Green's function treatment of the semi-infinite leads[85]. In summary, it still remains an 
interesting and unsettled question, whether intrinsic spin Hall transport can be experimentally 
observed in mesoscopic systems as studied in the above references. 

3.3 Spin Hall transport of holes 

We now analyze intrinsic spin Hall transport in p-doped III-V semiconductors and start with 
the case of a three-dimensional bulk system pioneered by Murakami, Nagaosa, and Zhang [12]. 

3.3.1 Three-dimensional bulk case 

We consider valence-band heavy and light holes governed by the Hamiltonian (12) and the 
conventionally defined spin current operator as described before. In this case a nonzero spin 
conductivity occurs if the direction of the spin current, its spin polarization, and the driving 
electric field are mutually orthogonal. For definiteness, let us assume the spin polarization to 
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point alone the z-axis with the electric field being in the x'y-plane. From the Kubo formula 
(38), the zero-frequency spin Hall conductivity can be evaluated as [59, 52, 86] 



6 71 + 272 



47r2 



72 
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f ^ dk- 



+ 



(45) 



where we have again assumed an infinite system at zero temperature, and 




71 T 272 



(46) 



are the Fermi wave numbers for heavy and light holes, respectively. Similarly to the approach 
leading to Eq. (41) for conduction-band electrons, disorder effects are taken into account via an 
effective relaxation time r. In the case here this approach is justified because vertex corrections 
can be shown to be absent for scattering potentials described by delta- functions [87]. Thus, 
for this type of disorder, the result (45) is exact in lowest order Born approximation [58]. 
The absence of vertex corrections in this case crucially relies on the fact that the underlying 
Hamiltonian (12) is not linear but of second order in the components of the momentum [87, 88]. 
The case of impurity potentials of longer spatial range was investigated in Ref. [89]. 

The remaining integral in Eq. (45) is elementary leading to a rather tedious expression 
which shall not be given here. However, we see that the value of the above integral is governed 
by the ratio of energy scale of the impurity scattering U/t and the "spin-orbit energy" 



n'^2{k'}f/m = 2£p2hi, 



(47) 



smce 



k% = ^2m£fhih'^ 



(48) 



is a typical wave number in the integration interval [59] . If ^/r 3> £so the spin Hall conductivity 
vanishes as 



TT 



+ 



n/Tj 71 



- 



2 / \ 2^ 
so \ {l2\ 



%/t J ' \h/T J \7l 



(49) 



where we have also assumed that the ratio 72/71 is small as it is usually the case [44]. In the 
opposite case h/r <^ Sgo one finds 



^f^(O) 




12 




(50) 



Here the the contribution in leading order is the result obtained in Refs. [52, 86] for a disorder- 
free system (up to some definitorial prefactor [59]). The expression given originally in Ref. [12], 
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however, differs somewhat from the above one due some approximation employed there [12]. 
Rcf. [90] contains calculations of the spin Hall conductivity in band structure models more 
general than Eq. (12). Numerical results based on ab initio band structure calculations were 
presented in Ref. [91]. A further numerical study of spin Hall transport within the Hamiltonian 
(12) in the presence of disorder was performed in Ref. [92]. 

In summary, spin Hall transport in p-doped bulk III-V semicondcutors is robust against 
disorder of not too large strength, but naturally breaks down if impurity effects are over- 
whelming the spin-orbit coupling [59]. 



3.3.2 Heavy holes in a quantum well 

Let us now turn to the case of spin Hall transport of heavy holes in p-doped quantum wells 
which was studied first by Schliemann and Loss [14]. An experimental observation of spin 
Hall effect in such a system was recently reported by Wunderlich et al. [15]. 
We consider the Hamiltonian (15) and a spin current 



p 3h 

Jz = TT^^" 

m I 



(51) 



for heavy holes with spin ±3/2 polarized along the growth direction of the well. Proceeding 
as above, one finds for the zero-frequency spin Hall conductivity [14] 
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(52) 



where r is again the momentum relaxation time. Similar to the three-dimensional bulk case, 
vertex corrections to the spin-Hall conductivity turn out to be zero for delta-function shaped 
scatterers [88], justifying the above approach. This result is again due to the fact that the 
spin-orbit coupling is not linear but of higher order in the particle momentum. 

The Fermi wave numbers entering Eq. (52) refer to the two dispersion branches (16) 
and can be expressed in terms of the particle density 



(53) 



as [14] 



kf - 



1 



2 2ma 



+ . 



/ 2ma\ 
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(54) 



Moreover, the longitudinal spin conductivities crf^^, ay^^ turn out out be identically zero. 

The remaining integral in the above expression (52) is elementary; however, it leads to 
a rather cumbersome expression which shall again not be given here. Analogously to the 
previous case, the energy scale of impurity scattering fi/r has to be compared with the "Rashba 
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energy" er = a(A;°)^, where k'j = y2me//^^ is the Fermi wave number for vanishing spin- 
orbit couphng, which is a typical value for k in the integration in Eq. (52). If the impurity 
scattering dominates over the Rashba coupling, h/r 2> er, the spin Hall conductivity vanishes 
with the leading order correction given by 

^a'('«-^|j!-^(('/-)^-(^-/)V«((|7)> (^^) 

where the Fermi wave numbers are given by Eq. (54). In the opposite case in S> Ii/t, the 
leading contribution to the spin Hall conductivity reads 



Note that this result for the spin Hall conductivity depends only on the length scale ma/h'^ 
of the Rashba coupling and the total hole density n, but not separately on quantities like 
the Fermi energy and the effective heavy hole mass. If ma/fi^ is small against the inverse 
square root of the total hole density (but still fulfilling er 3> ^i/r), the spin Hall conductivity 
approaches a value of 

4y' = 9^ • (57) 

This is the case if h/r <^ er <^ Ef. This above value should be compared with the universal 

value of e/S-TT found originally in Ref.[13] for electrons in a fully clean asymmetric quantum 
well. In this sense the hole spin Hall conductivity is enhanced by a factor of 9 compared to 
the naive result for electrons, which is partially due to the larger angular momentum of the 
heavy holes. 

Zarea and Ulloa studied the above system in the clean limit but with a perpendicular 
homogeneous magnetic field coupling to the orbital degrees of freedom of the holes but not 
to their spin [93], an investigation analogous to the one by Rashba on n-doped quantum 
wells already mentioned [71]. Again it is found that, for an infinite system, the case of zero 
magnetic field and the limit of vanishing magnetic field do not coincide [93]. The details of 
this effect, however, seem to be somewhat different from the observations made in Ref. [71] 
and need further study. In any case, influence of a magnetic field coupling to the orbital 
degrees of freedom only should only be appreciable if the field is strong enough to produce 
typical cyclotron radii being of order of the system size or smaller. Therefore, arbitrarily small 
fields cannot be expected to have an effect in real experiments. Another study of spin Hall 
transport in the presence of a perpendicular magnetic field was performed in Ref. [94] where 
a more general band structure Hamiltonian was used [88] 

As seen above, spin Hall transport of heavy holes in a quantum well is robust against 
disorder effects, differently from the situation for electrons in n-doped wells. This effect is 
due to the different functional form of the effective spin-orbit coupling and was also confirmed 
numerically by Nomura et al. who performed a careful comparison between these two systems 
[75]. In a subsequent study, the edge- spin accumulation caused by the spin current was 
investigated numerically [95] . Further investigations of disorder effects can be found in Ref. [96] 
where an approach based on nonequilibrium Green's functions was used. 

Moreover, several groups have studied numerically tight-binding models for two-dimensional 
hole systems coupled to semi-infinite leads. [81, 97, 98]. These investigations are analogous 
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in spirit and tcclinical approach to the numerical work on models for n-doped systems men- 
tioned earlier [76, 77, 78, 79, 80, 82, 83, 84]. In particular, the numerical results on heavy-hole 
systems [81, 98] in the limit of low disorder confirm quantitatively the the enhanced spin con- 
ductivity (57) obtained analytically in Ref. [14]. Finally, many of the abovementioned mainly 
numerical studies on p-doped quantum wells were inspired by the experiments by Wunderlich 
et al. which we will discuss in section 4. 

3.4 Spin Hall effect in other systems 

Intrinsic mechanisms of spin Hall transport in n-doped bulk III-V semiconductors where inves- 
tigated by Bernevig and Zhang [99, 100]. Here the leading contribution to spin-orbit coupling 
is given by the bulk Dresselhaus term (3) being of third order in the electron momentum. On 
the other hand, Engel, Halperin, and Rashba have studied extrinsic spin Hall effect in such 
systems [101]. We will discuss these issues in section 4 in the context of the experimental 
results by Kato et al. [16]. 

Spin Hall effect in graphene, i.e. single planes of graphite, was investigated by Kane and 
Mele [102]. Using a theoretical picture similar to the edge-state theory of the charge quantum 
Hall effect, these authors propose a spin current at the edges of a graphene sheet. Work 
following these theoretical predictions include Refs. [103, 104]. Finally, Shchelushkin and 
Brataas considered spin Hall transport in normal metals due to extrinsic mechanisms [105]. 

4 Detection of spin Hall transport: Experiments and Propos- 
als 

We now turn to experimental investigations of spin Hall effect in semiconductors. The studies 
already carried out and many of the experimental proposals found in the literature use the spin 
accumulation caused by the spin current for detecting spin Hall transport. Using a simple spin 
diffusion model, this spin accumulation is expected to decay towards the bulk of the sample 
on a length scale given by the spin diffusion length [106]. The latter quantity is determined 
by the semiconductor material, but possibly also by further details of the sample and the 
experimental setup [107]. Further recent theoretical studies on spin accumulation caused by 
the spin Hall effect include [95, 108, 109, 110, 111, 112]. 

Kato et al. have studied spin Hall transport in n-doped bulk epilayers of GaAs and InGaAs 
with a thickness of 2//m and 500nm, respectively [16]. The electron density in both samples 
was n = 3 • lO^^cm"^. The spin Hall effect was detected via the optical technique of Kerr 
rotation microscopy in the presence of an external magnetic field in a Hanle-type setup. The 
Kerr rotation signal as a function of the applied magnetic field could be fitted well by a 
Lorentzian where the spin lifetime Tg entered as a fit parameter. By scanning over the sample, 
the authors obtained spatial profiles of the spin density along the direction perpendicular to 
the applied electric field. Fitting this data to solutions of the spin diffusion equation, the 
spin diffusion length was extracted. The values for this quantity lie between two and four 
micrometers and are, within the error bars, insensitive to the electric field varying between 
zero and 25mV;um~^. Combining this data with results for the spin lifetime, the authors 
inferred a spin conductivity of about 0.5r2~^m~^, where the latter quantity was converted to 
units of charge transport by multiplying with a factor of e/h. 
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The samples investigated by Kato et al. arc n-dopcd and clearly in the bulk regime. There- 
fore, the intrinsic spin-orbit coupling to conduction-band electrons is dominated by the bulk 
Dresselhaus term (3) which was studied by Bernevig and Zhang as a possible intrinsic mech- 
anism for spin Hall transport. [100]. The authors apply their results to the experiments by 
Kato et al. [16]. However, the agreement between theoretical predictions and the experimen- 
tal findings is certainly not very convincing. Moreover, Kato et al. find only a very negligible 
dependence of the spin Hall transport on strain applied to the system. This observation also 
strongly disfavors an intrinsic mechanism. In fact, Engel, Halperin, and Rashba [101] have 
developed a theory of extrinsic spin Hall transport in GaAs based on impurity scattering and 
found reasonable agreement with the results by Kato et al. [16]. 

Further results on this type of experiments were reported by Sih et al. from the same 
research group for the case of n-doped GaAs quantum wells [113] (as opposed to the bulk 
epilayers discussed so far) . The quantum well here has a width of 75 A with a sheet density of 
n = 1.9- lO^^cm"^ and a mobility of = 940cm^/Vs. From the analysis of their Kerr rotation 
data Sih et al. conclude that the signatures of spin Hall transport seen in their experiment 
are also most likely due to an extrinsic mechanism. 

We now turn to spin Hall transport of holes. Wunderlichei al. have investigated spin Hall 
effect in a p-doped triangular quantum well which is part of a p-n junction light emitting 
diode [15]. The quantum well has a sheet hole density of ri = 2 • lO^^cm"^ which is for the 
still low enough, for the given sample geometry, such that only the first heavy hole subband is 
occupied. Thus, the intrinsic spin-orbit coupling to these heavy holes can be expected to be 
governed by the Hamiltonian (15), leading to a spin Hall conductivity given by Eq. (52). The 
spin accumulation at the edges of the well is detected by the circular polarization of the light 
emitted from the diode. In a subsequent publication, the authors presented further details on 
this technique [114]. In summary, Wunderlich at al. conclude that the spin Hall transport 
seen in their results is most likely of intrinsic nature, i.e. due to Rashba spin-orbit coupling as 
described by the Eq. (15), although further work is needed to fully establish this conclusion 
[15]. 

Let us now mention some proposals for further possible experimental investigations of spin 
Hall transport. Hankiewicz et al. have considered an H-shape two-dimensional electron system 
using the Landauer-Biittiker formalism [76]. The sample consists of two, say, horizontal bars 
connected in their centers by a shorter vertical bar. A voltage applied along one of the long 
bars should drive a spin current through the shorter connecting bar into the other horizontal 
bar. There this spin current should manifest itself by a voltage along the bar which can 
be calculated by inverting the spin conductivity tensor. Numerical calculations support the 
feasibility of this experimental approach [76]. 

A scheme to determine the spin Hall conductance purely via measurements of charge trans- 
port was put forward by Erlingsson and Loss [115]. The authors study a planar four-terminal 
setup. Using conventional scattering formalism, they express the spin Hall conductance in 
terms of voltages, charge conductances, and charge current noise quantities. This allows, in 
principle, to infer the spin Hall conductance just from electric measurements, avoiding the 
necessity of any magnetic or optical element [115]. 
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5 Conclusions and outlook 



We have given an overview on recent theoretical and experimental developments concerning 
spin Hall transport in semiconductors. This phenomenon has certainly been over the last 
years one of the most intensively worked on topics in the solid-state community, and research 
efforts still continue to grow. 

The theoretical situation regarding spin Hall effect in the two-dimensional electron gas 
with spin-orbit coupling linear in the momentum is by now well settled: There is no spin 
Hall effect in an infinite system if any kind of dissipation mechanism is present. This is, 
however, not a statement about finite mcsoscopic systems, and further both theoretical and 
experimental work is to be expected here. On the other hand, p-doped quantum wells, i.e. 
the two-dimensional hole gas, appears to be a particularly attractive system. 

A challenge for future theoretical work is certainly to extend many present results to the 
case of finite temperature, and, more importantly, include the Coulomb interaction between 
charge carriers. The last point is addressed by only very few papers so far [116, 117]. A 
further recent theoretical development are studies on the zitterbewegung of electron and hole 
wave packets due to spin-orbit interaction in semiconductors [118, 119]. These systems offer 
the possibility to experimentally detect the relativistic phenomenon of zitterbewegung which 
appears to be quite inaccessible in the case of free electrons [118, 119]. 

Future challenges for experiments on spin Hall transport include the clarification and 
discrimination of extrinsic and intrinsic mechanisms. Moreover, from both a theoretical and 
an experimental point of view, it is undoubtedly desirable to develop setups and techniques 
which allow to detect spin Hall transport independently from spin accumulation. 
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